Individu: Poisson
Population: Le lez
Échantillon: 100 poissons dans l’aquarium
Taille: 100
Variable: Couleur du poisson
Ici la variable a bien les modalités suivantes:
Incompatible: Un poisson ne peut avoir qu’une seule couleur dans cette étude pour simpiflier la representation
Exhaustive: Visiblement le lez ne contient que des poissons rouges, vert ou bleus, dans les faits je ne suis pas certains que ce soit le cas
Sans ambigité: En partant du principe que les 3 couleurs de poissons sont rouge vert et bleu, il n’y a pas d’ambiguité sur la couleur du poisson, dans les faits on pourrait prendre la couleur la plus presente sur le poisson
set.seed(2002)
lez <- c(rep("red", 2000),
rep("green", 3000),
rep("blue", 5000))
n_exp <- 100 ## Nombre de fois où je tire un poisson
echantillon <- sample(lez, ## échantillonne les poissons
n_exp, ## n_exp fois
replace = FALSE) ## sans remise
Résumé
res <- table(echantillon)
res
## echantillon
## blue green red
## 52 22 26
La fréquence
freq <- prop.table(table(echantillon))
freq
## echantillon
## blue green red
## 0.52 0.22 0.26
On a eu une proportion de poisson bleu plus elevés, ce qui semble logique vu qu’il y a plus de poisson bleus dans le lez, notre “chance” de piocher un poisson bleu est donc plus importante que pour un poisson vert ou rouge.
Graphique
barplot(res, col = names(res) )
set.seed pour
produire un document reproductible.Quand on compile plusieurs fois on obtient toujours les memes résultats car la seed est la meme
set.seed(12)
n_exp <- 10 ## Nombre de fois où je tire un poisson
echantillon <- sample(lez, ## échantillonne les poissons
n_exp, ## n_exp fois
replace = FALSE) ## sans remise
res_10 <- table(echantillon)
barplot(res_10, col = names(res_10) )
Ici on remarque que le graphique ne represente pas bien du tout notre population, c’est du au fait que l’échantillon soit trop petit
n_exp <- 100 ## Nombre de fois où je tire un poisson
echantillon <- sample(lez, ## échantillonne les poissons
n_exp, ## n_exp fois
replace = FALSE) ## sans remise
res_100 <- table(echantillon)
barplot(res_100, col = names(res_100) )
Ici on remarque que le résultat represente déjà bien mieux la population pour les poissons bleus, pour avoir une idée plus précise des poissons verts et rouge il faudrait un échantillon encore plus grand